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Abstract. Wc address binary optical communication in single- mode and entangled 

I quantum noisy channels. For single-mode we present a systematic comparison between 

' direct photodetection and homodyne detection in realistic conditions, i.e. taking into 

^ , account the noise that occurs both during the propagation and the detection of the 

' signals. We then consider entangled channels based on twin-beam state of radiation, 

■ and show that with realistic heterodyne detection the error probability at fixed channel 
On ! energy is reduced in comparison to the single-mode cases for a large range of values of 

' quantum efBciency and noise parameters. 

o : 

■ 

, 1. Introduction 

^ : 

^ ' Classical information may be conveyed to a receiver through quantum channels. To this 
O^' aim a transmitter prepares a quantum state drawn from a collection of known states and 
^ ■ sends it through a given quantum channel. The receiver retrieves the information by 
^ ■ measuring the channel, such to discriminate among the set of possible preparations, and 

■ to determine the transmitted signal. The encoding states are generally not orthogonal 
and also when orthogonal signals are transmitted, they usually lose orthogonality 
because of noisy propagation through the communication channel. Therefore, in general, 
no measurement allows the receiver to distinguish perfectly between the signals ^ |2] 
and the need of optimizing the detection strategy unavoidably arises. 

In binary communication based on optical signals, information is encoded into two 
quantum states of light. Amplitude modulation-keyed signals (AMK), consist in two 
states of a single-mode radiation field, which are given by Qj = j = 1,2, with 

li^i) = \i^o) and \ip2) = D{a)\ipo), where \ipo) is a given seed state, usually taken as the 
vacuum, D{a) = exp(a;a"'' — a*a) denotes the displacement operator, and the complex 
amplitude a may be taken as real without loss of generality. The reason to choose 
AMK signals lies in the fact that displacing a state is a simple operation, which is 
experimentally achievable by a linear coupler and strong reference beam Another 
binary encoding, based on displacement, is given by phase shift-keyed signals (PSK) 
\4>i) = D{—a)\ipo) and \ip2) = D{a)\ipo)- In the following, we will refer to the AMK 
situation only, though all the results hold also for PSK. The price to pay for such a 
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convenient encoding stage is that, for any choice of the seed state, the two signals 
are always not orthogonal, and thus a nonzero probability of error appears in their 
discrimination |3]. If we consider equal a priori probabilities for the two signals, the 
optimal quantum measurement to discriminate them with minimum error probability is 
the projection- valued measure {Mj}j=i^2, Mi + M2 = I, corresponding to P 

M, = ^T[(-yA,]|A,)(A,|, (1) 

k 

where \Xk) is an eigenstate of the hermitian operator A = p2 — pi with eigenvalue 
Afc, and T[x] is the unit step function, which is zero for negative x, one for positive 
and T(0) = ^. The probability of inferring the symbol j when i is transmitted is 
given by P{j\i) = Tr{f)j Mj}, j,i = 1,2, such that the average error probability in 
binary communication is given by = | {P(l|2) + P(2|l)}. The minimum of Pe, 
corresponding to the optimal measurement (^, reads as follows 

i'e = ^ (1 - VT^VMhW) , (2) 

and is known as the Helstrom bound 0. In particular, for a pair of AMK signals we 
have 



^^^Wizi5^, (3) 

where with we denote the average number of photons in the channel per use, i.e. 
N = |Tr[a"''a(po + Pa)] = + '^o, where uq = {ipo\a^ a\ipo) is the average photon 

number of the seed state. For the sake of brevity, we will refer to also as to the 
energy of the channel. Notice that for a pair of PSK signals the same bound in equation 
^ holds; however, the expression for N is now given by = |ap + rio- 

Binary communication has been the subject of much attention, mostly concerning 
the design and the implementation of optimal quantum detection processes, to 
distinguish nonorthogonal signals with reduced error probability, possibly approaching 
the Helstrom bound given in equation The relevant parameter in this optimization 
is the energy of the channel, which itself limits the communication rate of the channel. 
After the pioneering work of Helstrom pp, a near-optimum receiver for AMK signals 
based on direct detection was proposed in j3] , whereas an optimum receiver approaching 
the minimum error probability (based on photon counting and feedback) has been 
suggested in jH]. More recently, various efforts has been made to find out optimum 
detection operators and decision processes for more general signals and in presence of 
noise [3 IHl El HD] • Indeed, when one has at disposal a given set of quantum signals, 
the problem becomes that of finding the optimal receivers E] and detection schemes 
[TT| IT^ and to compare their performances with those of realistic detectors. Following 
this way, some studies were made on the effects of thermal noise on the optimum 
detection for a coherent AMK channel [T^ . 

In this paper we focus our attention on protocols for binary communications where 
both AMK signals and receivers can be realized with current technology. As we will 
see, the various sources of loss and noise can be described as an overall Gaussian noise. 
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Our analysis allows to unravel the different contributions and to compare receivers in 
realistic working regimes. The purpose is twofold: on one hand we perform a systematic 
comparison between direct and homodyne detection in presence of noise during the 
propagation and the detection stages, in order to find in which working regimes a 
receiver should be preferred. On the other hand, we show that binary communication 
can be improved by using achievable sources of entanglement and realistic heterodyne 
receivers. Indeed, it has been recently shown that in ideal conditions (perfect detection 
and noiseless propagation) entanglement improves the performances of a binary channel, 

1. e. it reduces the error probability in the discrimination of the symbols ^1 IT^ . 
Motivated by these results, we investigate the error probability of entangled channels 
in realistic conditions, taking into account the unavoidably noise that occurs during the 
propagation and the detection. Since we are interested in assessing entanglement as 
an effective resource, we compare entangled channels with the corresponding realistic 
single-mode channels. 

The paper is structured as follows. In Section [21 we address single-mode channels 
that uses direct or homodyne detection as receivers, and compare the corresponding 
error probabilities both in ideal and realistic situations, i.e. in presence of noise. In 
section El we describe a binary communication scheme based on entangled twin-beam 
state of radiation that employs multiport homodyne or heterodyne detection in the 
measurements stage. As we will see, there are regimes where the error probability is less 
than in a single mode channel, also when the noise affects propagation and detection. 
Finally, in section |31 we summarize our results giving some concluding remarks. 

2. Single Mode Communication 

2.1. Direct detection - Ideal case 

A scheme based on direct detection, to discriminate the set {|0), can be 

implemented as in figure ^ ISj- It consists of a beam splitter (BS), in which the state 
to be processed, either |0) or \a), is mixed with a given coherent reference state, say 
1/3). The outgoing mode is subsequently revealed by on/off photodetection, i.e. by a 
detector which checks the presence or absence or any number of photons. 

The operator describing the action of a BS on the modes a and b of the field is 



where r = cos^ is the BS transmissivity and a, and b, 6^ are the annihilation and 
creation operators for the two modes, respectively. If the input state is \a)a\l3)b = 
\ci)a'^ \l3)b, |«)a and \P)b being coherent states, the output state is given by \c()a\P)b — »■ 
I a cos (j) + j3 sin (j)) a \ — a sin -|- cos (j))b. By choosing the amplitude of the reference as 
P = —a/ tan0, we obtain, at the output, the vacuum when the input state is \a), and 
I — Q!cos0) for vacuum input, in formula 





a) |0) 



and 




(5) 
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In order to discriminate the two input signals, one performs a simple ON / OFF 
photodetection: when the output is the vacuum the detector doesn't click, otherwise 
it clicks. This measurement is described by the probability operator-valued measure 
(POVM) {Ho, Hi}, where Ho = |0)(0| and Ho + Hi = I, i.e. we assumed unit detector 
efficiency. The error probability, K^, is defined as: 

K, = ^{K{0\a) + K{am, (6) 

where K{0\a) and K{a\0) are the probabilities of inferring that the input state is |0) 
when it is actually |a) and vice versa. In our case 

KiO\a) = Tt{U^ p„ ® pp ul Hi ® 1} , (7) 

K{a\0) = Ti{U^ Po ® p/3 ul Ho ® 1} , (8) 

where p^ = p. = 0,a,p. We obtain 

K(0|a) = |(0|ni|0)p = (9) 

K{a\0) = |(— acos0|no| — acos</))|^ = exp{ — cos^ </>}, (10) 
such that Eq. (P) rewrites as 

exp(— |ap cos^ 0) exp(— 2A^ cos^ 0) , , 

Ke = = . (11) 

with = TT{a^a{po + pa)} = Notice that in the limit |ap ^ 1 (relevant 

for classical communication) and cos^ ^ 1, we have 2Pc- This is usually 

summarized by saying that the measurement is asymptotically near optimal. 



2.2. Direct detection - Noise in propagation and detection 

In this section we take into account the effects due to the noise that occurs in the 
propagation and the detection of the signals. We model the propagation in a noisy 
channel as the interaction of the single-mode carrying the information with a thermal 
bath of oscillators at temperature T. The dynamics is described by the Master equation 

^ = {r(l + M)L[a] + TML[a^]}p, , (12) 

where pt = pit) is the density matrix of the system at the time t, T is the damping 
rate, M = [q^/^bT _ jg ^.j^g number of thermal photons with frequency uj at 
the temperature T, and L[0] is the Lindblad superoperator, L[0]pt = OptO"^ — 
^O'^Opt — \ptO'^0. The term proportional to L[a] describes the losses, whereas the 
term proportional to L[a'^] describes a linear phase- insensitive amplification process. 
In other words, we are taking into account the unavoidable dissipation and in-band 
amplifier noise. We are not considering other sources of noise such as cross-talk and 
inter-symbol interference. The Master equation p2|) can be transformed into a Fokker- 
Planck equation for the Wigner function 

W,{0 ^ W[pM) = -2 e^^*"'^*' Tr {p, D{\)} = - exp {-2|C - p\'] , (13) 

71 J n 
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with ( & C, = /i = 0,a, and -D(A) is the displacement operator. Using 

the differential representation of the Lindblad superoperator ^HJ^], the Fokker-Planck 
equation associate to equation (fT^ is 



>,y) (14) 



where we put ( = x + iy and (C) is the Wigner function of the system at time t, for 
initial state p^. The solution of equation ()14|) can be written as the convolution 



W,,t{x, y) = J J dx'dy' W,,o{x', y') Gt{x\x') Gt{y\y') , (15) 
where the Green function Gt{xj\x'j) is given by 

1 f (x, -x;.e-^r*)2] 

with = |(M + i)(l — e~^*). Using the equations (jl5p and (jl6p . we arrive at 

^..*(C) = exp |-KZ|^Z!!^ I , (17) 

with A|^p = |[1 + 2M(1 — e~^*)]. The Wigner function in equation (fTTj) corresponds to 
the density matrix of a displaced thermal state (THj 

p^{t) = pM' = D{a') VM' D\a') , (18) 

where a' = a e~^^* and vm' is a thermal state 

1 / M' y'" , , 

^M' = --TTT --tt: (19) 



M' \l + M\ 

with M' = M(l — e~'"*) average number photons. 

Equation (|THjl describes the signal arriving at the receiver (figure [T]). Let us now 
consider the noise in the detection stage. At first we have to choose the reference state. 
Since this receiver is based on the interference between the signal and the reference, it 
turns out that, in presence of propagation noise, the optimal reference is the coherent 
state 1/9'), j3' = —a'/ tan0. Moreover, if on/off detection is not ideal, we must consider 
the finite detector efficiency rj. In this case, the POVM describing the measurement is 

Ho (77) + Hi (77) = I, where 

00 

Uo{v)= Y,{l-vr\n){n\. (20) 

n.=0 

In order to evaluate the detection probabilities we use the fact that the trace between 
two operators, Oi and O2, can be written as the phase-space integral ^\ 

Tr{Oi02} = 7r J d\W[0,]iC)W[02]iC) , (21) 
where the Wigner function of a generic operator O is defined as 

W[0]{C) = 1^ j d'Ae'^^*^^*^Tr{OZ^(A)}. (22) 
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The Wigner of the POVM element Ho is then 

with A2 = (1 -r/)/(2r/). 

The expression (jSI) for probabihty A'^^r,A/(«|0) to infer \a) when |0) is sent is 
modified as follows 

i^r,,r,M(a|0) = TiiU^UM' ® (/3'|?7j no(ry) ® 1} . (24) 
By equation pT|l we have 

i^,,r,M(a|0) = tt' jd'-fd'XW[U^UM' ® t/i](7, A) iy[no(r/)](7) , (25) 

which, using 

W[U^ UM' ® 1/3') [/J] (7, A) = W[uM'] (7 cos + A sin 0) 

X PF[|/?')(/3'|](-7sin0 + Acos0), (26) 

leads to 

A-,..M(aio) ^ . (27) 

1 + ?7M(1 — e ^*)cos^0 
On the other hand, for the probability i^'^,r,M(0|a) we have 

i^,,r,M(0|a) = TT{U^D{a')uM'D\a') ® |/5')(/3'|?7i ni(r/) ® 1} , 

= l-TT{U^D{a')uM' D\a') ® no(r/) ® 1} , 

?7M(1 - e"^*) cos^ 
~ 1 + r/M(l -e-r*)cos2 ' 

and, finally, the overall error probability in presence of noise is 



28 



/.r/V FMl .M(l-e-^OcosV + exp 

^' ^' = 2[l+r/Af(l-e^r.),,,2^] ' (29) 

which reduces to equation (jllj) in the limits 77 ^ 1 and F, M — > 0. Notice that the 
result of equation corresponds to the presence of an overall Gaussian noise with 
parameter cj^ = riM{l — e~'"*) plus dissipation of the signal {a ae~2^*). Our analysis 
allows to unravel the different contributions to a. 



2.3. Homodyne detection - Ideal case 

An alternative receiver for single mode communication is provided by homodyne 
detection, which offers the advantage of amplification from a local oscillator, avoiding 
the need of single-photon avalanche photodetectors [SO]- A schematic diagram of the 
balanced homodyne detection is shown in figure El here the signal interferes with a local 
oscillator (LO), i.e. a highly excited coherent state, in a balanced BS [this corresponds to 
put = 7r/4 in equation (jl))]. After the BS the two modes are detected and the difference 
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photocurrent is electronically formed. For unit quantum efficiency of photodiodes, the 
POVM of the detector is = |a;)(x|, with 

= - e > — } \n) (30) 

^ ' n=0 

being an eigenstate of the quadrature operator x = |(a + a''') of the measured mode. In 
equation (|3Up Hn{x) denotes the n-th Hermite polynomials. The probability density 
p{x\fi) of obtaining the outcome x from homodyne detection with input state 
/i = 0, a, is 



p{x\fi) = Tt{\i2) n J = \l=- exp{-2(x - fif} . (31) 

In equation (jHT|) a is assumed as real. Equivalently, if a G C, the same result may be 
obtained by measuring a suitable quadrature x^, = |(a^'e'^ + ae'"^), with ip = arg(a). 
The minimum error probability for the homodyne receiver is given by 

H, = ^{H{0\a) + H{a\0)} 



1 1 /"'^/^ 

dxp{x\a) + / dxp{x\0) 

-oo Ja/2 



Erf 



(32) 



2 

_ 1 
~ 2 

where H{0\a) and H{a\0) are the probabilities of inferring the signal |0) when it is 
actually |a) and vice versa. Erf [a] = dCe~^^ denotes the error function. Notice 

that, in general, the error probability depends on the choice of a threshold parameter 
A, i.e. 

"A 



^^o(A) = I 



/A p+oo -\ 

dxp{x\a) + / dxp{x\0) > 
■oo J A ) 



in 



i (Erf |y2Al + Ed[-/2(a - A)]) } ■ (33) 

In our case this probability is minimized when A = a/2, thus leading to the result in 
equation (j22I). In the limit ^ 1, equation (j22)) reduces to 

e"^ 

He ~ ^= . (34) 

Homodyne detection provides a better discrimination of the signals than direct 
detection, i.e. < if the energy of the channel is below a threshold that 

monotonically increases as the transmissivity of the BS in the receiver decreases. In 
the limit of r — > 1 in direct detection, we have He < ioi N < 0.77 while, as an 
example, for r = 0.9 we have He < Ke for < 1.10. 

2.4- Homodyne detection - Noise in propagation and detection 

As already shown in section YI7I[ if noise affects the propagation of the signal, the state 
arriving at the homodyne receiver is no longer a pure state, and is given in equation 
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(I18p . Moreover, when also homodyne detection is not ideal, the POVM of the receiver 
is a Gaussian convolution of the ideal POVM 

where = (1 — ?7)/(4?7), and r] is the quantum efficiency of both photodiodes involved 
in homodyning (we assume that they have the same quantum efficiency). The Wigner 
function of na;('7) is given by 

1 r (Re[C]-x)^ ) 



Taking into account all the sources of noise, the probability density of equation (j3T| 
becomes 

P,,r,M(a:|/i) = Tr{D(/i') um' D\^i') liM) , (37) 

with /i' = /i e~^^*, /i = 0,a, and um' given by equation (fTUI). In this way, using the 
Wigner functions and thanks to equation ()21|) . the error probability reads as follows 



i7e(iV,r/,r,M) = ^<|l-Erf 



/r/iV e 



V^l + 2r/M(l - e-r*) 
which, in the limit r/A^ e~^* ^ 1 + 2?7Af (1 — e"'"*), reduces to 



(3^ 



. '1 + 2vM{l - e-r*) exp {- | 



In the next section we compare with the corresponding error probability in direct 
detection. 



2.5. Direct vs homodyne detection 

In order to individuate the working regimes where homodyne detection provides better 
performances, i.e. lower error probability than direct detection, we define the following 
quantity 

Ae{N, r/ken, Vhom, T, M = 1 - — — — — , (40) 

Kc{N, r/ken, r, M) 

where ?7kcn and ?7hom are the ON / off and homodyne detection efficiencies, respectively. 
When Aq > homodyne receiver's error probability is the lowest. In figureElwe plot 
as a function of the energy of the channel N for different values of the other parameters. 
For given values of r/ken; Vhom, T and M we have two different thresholds for the energy 
channel, namely Nl^^°^\ri]^en,Vhom,^, M), j = 1,2, such that Ac = 0. For < A^th,°™'* 
we have a small interval where homodyne detection should be preferred to direct one; 
as increases, we find a window {N^^"^^ < N < N^^"^^), where > and, finally, 
a last region for > N^^^^ where homodyne detection returns to be definitely better 
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(see table 121). This result is due to the presence of thermal noise {i.e. T ^ 0). In fact, 
as one can easily see from equations ()29|) and (j^HI) one has, independently on M 

hm K^[N, T^kcn, r, M = -y— — — — , 41 

A^-^oo 2[1 + r^kcn^ll — e ^*)cos^0J 

lim H,{N, r/hom, T, M) = . (42) 

In summary, homodyne detection provides better results for either small or large 
values of the channel energy A^, whereas for intermediate values of the optimal choice 
is represented by direct detection. The width of this intermediate region decreases as 
the noise increases, i.e. as the value of both F and M increases. We therefore conclude 
that homodyne detection is a more robust receiver in presence of noise. As concern 
quantum efficiency, we have, as one may expect, that the performances of each detector 
improves increasing the corresponding rj. 



3. Binary communication in entangled channels 

Entanglement is a key feature of quantum mechanics. The quantum nonlocality due to 
entanglement has been, in the last decade, harnessed for practical use in the quantum 
information technology [211 122] • Entanglement has become an essential resource for 
quantum computing j22j, quantum teleportation [23], dense coding [TH], and secure 
cryptographic protocols |22] as well as for improving optical resolution |21j, spectroscopy 
[23], and general quantum measurements [2^]- Here we show how entanglement, and in 
particular entangled states that can be realized by current optical technology, can be 
used to improve binary communications, i.e. to reduce the error probability at fixed 
energy of the channel. 



3.1. Heterodyne detection - Ideal case 

Binary optical communication assisted by entanglement may be implemented using 
twin-beam (TWB) state of two modes of radiation ^3]. Schematic diagrams of some 
possible implementations are given in figure |3] In the Fock basis the TWB writes as 
follows 

|A)) = v/r^5^ A'M^)|n), (43) 

n 

where |A| < 1 and, without loss of generality, it may be taken as real (A is sometimes 
referred to as the TWB parameter). TWB is the maximally entangled state (for a 
given, finite, value of energy) of two modes of radiation. It can be produced either 
by mixing two single-mode squeezed vacuum (with orthogonal squeezing phases) in a 
balanced beam splitter [2S| or, from the vacuum, by spontaneous downconversion in a 
nondegenerate optical parametric amplifier (NOPA) made either by type I or type II 
second order nonlinear crystal [2Z]. Referring to the amplification case, the evolution 
operator reads as Ur = exp {r (a^6^ — a6)} where the "gain" r is proportional to the 
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interaction-time, the nonlinear susceptibility, and the pump intensity. We have A = 
tanh r , whereas the number of photons of TWB is given hjNx = 2sinhV = 2AV(1-A2). 

The two signals to be discriminated in a AMK encoding are given by \ipo)) = 1-^)) 
and I'lpa)) = [-Da(«) ® I]|A)), where the displacement operator Da{a) is acting on one 
of the modes, say a, of the TWB. The energy of the TWB channel, i.e. the average 
photon number per use, is given by = |Tr{ (a''"a C?) I + 1® 6^6)(po + Pa)} = A'a + ^|ap, 

with po = \i'o)){{i'o\ and pa = \^a)){{^a\- 

The error probability for the ideal discrimination between the two states \ipo)) and 
li^a)) reads as follows 



^ _ - v^l-|((^0|^a))P 

^ l-v/l-exp{-27V(i-/3)(i + /?iV)} ^ ^^^^ 

where (3 = Nx/N is the fraction of the channel energy that is used to establish the 
entanglement between the two modes. Probability ()44|) is minimum for (3 = {N — 1)/2N, 
when > 1, and for /3 = when < < 1. In summary one has 

Qe = Po {N< 1) (45) 



1- Jl-exp{-i(l + Ar)2} 
Qe = ^- {N>1), (46) 

where Pe, given in equation (j21), is the minimum error probability for single-mode AMK 
signals. Equations and say that Qc < Pc VA^, i.e. that the use of entanglement, 
at least in the ideal situation considered so far, never increases the error probability, 
and it is convenient if the photon number of the channel is larger than one. 

In order to see whether this result holds also in practice, it is necessary to find out 
a realistic receiver able to discriminate \ipo)) and \ipa)), and to discuss its performances 
in presence of noise. As concern detectors, we may use either multiport homodyne 
detection, if the two modes have the same frequencies [2H1; or heterodyne detection 
otherwise j2H]- Both these detection schemes allow the measurement of the real and 
the imaginary part of the complex operator Z = a — b"^ . In figured we have referred to 
to eight-port homodyne detection; however all the results also hold for other multiport 
homodyne schemes and for heterodyne detection. Each outcome from the measurement 
of Z is a complex number z and the POVM of receiver is given by 11^ = ^\z)){{z\, where 

\z)) = [Da{z) ® I] 5^ \n)\n) = [I ® D,{-z*)] \n)\n) . 

n n 

As already discussed for homodyne detection, we may take the amplitude a as real. 
In this case a suitable inference rule to infer the input state from Z-data involves the 
real part of the outcome as follows 

MA > A ^ l^a)) , 
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where A is a threshold value, which should be chosen such to minimize the probability 
of error 



Re = ^{R{0\a) + Ria\0)} 



(47) 



where R{0\a) and R{a\0) are the probabilities to detect \ipo)) when \ipa)) was sent and 
vice versa. The heterodyne distribution conditioned to a displacement D{a) = Da{a)^I 
is given by the probability density 



with fi 
Therefore we have 



r(.|/i) = |((.|D(/i)|A))p 

(1 - A)/(l + A) 



exp 



l/i 



A! 



0, a, and A\ 



{VW+2 - v^)/(v/iVA + 2 + 



R{0\a) 
R{a\0) 



dx 
dx 



dyr{z\a) 
dyr{z\0), 



where z = x + iy, and the error probability ()47|) becomes 



Re. 



l-i(Erf 



A 

a7 



+ Erf 



a 



A' 



A, 



(48) 
^). 

(49) 
(50) 

(51) 



-Ro in equation (jKT|) is minimized by choosing A = a/2, thus leading to 



Rp 



1-Erf 



1 a 
2AI 



At this point, the error probability Rc can be further minimized by tuning the 
entanglement fraction j3. By substituting the expression for the amplitude a = 
^2N{1-P) and the variance A 



(vWT2 - vW)/(\/^ 



-Re 




Erf 



'2N{l-/3)iVW- 



vm 



2 + A/5iV), we obtain 
(52) 



The optimal entanglement fraction, which minimizes -Re, is given by 



2(1 + N) • ^^^^ 

In figure El we report the resulting expression for the error probability, compared with 
the corresponding single-mode error probabilities and He- As the channel energy 
increases, the heterodyne error probability decreases more rapidly than the single-mode 
ones. For a channel energy < 0.79 homodyne detection gives the best performances, 
whereas as the channel energy increases the best results are obtained by direct detection 
(0.79 ^ N < 4.46) and heterodyne detection (A^ > 4.46). These results are summarized 
in tabled Notice that for A^ > 5.2 we have -Re < -Pc, i-G. TWB heterodyne channel 
provides better performance even than ideal single-mode channel. 
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3.2. Heterodyne detection - Noise in propagation and detection 

At first we consider the noise occurring during the propagation of the TWB or its 
displaced version. This is described as the coupling of each mode of the TWB with a 
thermal bath of oscillators at temperature T. The dynamics is then described by the 
two-mode Master equation 

^ = {r(l + M)L[a] + r(l + M)L[h] + TML[a^] + TMLp]] pt , (54) 

where pt is the density matrix of the bipartite system and the other parameters are as in 
equation (fT^. The terms proportional to L[a] and L\b] describe the losses, whereas the 
terms proportional to L[a^] and L\b)^] describe a linear phase-insensitive amplification 
process. Of course, the dissipative dynamics of the two modes are independent on each 
other. 

The Master equation can be reduced to a Fokker-Planck equation for the two- 
mode Wigner function W^{^, Q = W[pfj]{^, () of the system, 

W,{^,0^^, jd'x jd'Xe^^*-^'^e<''~^''TT{p^Da{x)<8)D,{X)} , (55) 

where ^, C G C, = I'ipfi)) {{i'fj.l, p = 0,a, and Dj, j = a, b, is the displacement operator 
acting on mode j. Using the differential representation of the superoperator in equation 
fl54|) . the corresponding Fokker-Planck equation reads as follows 

dtW,M^ C) = ^ |E (^-.^^- + ^y^y^) + + 5,%,) I w,M. C) , (56) 

I j=i i=i J 

where ^ = Xi + iyi, C = ^2 + W2- The solution of equation can be written as 

W,M^ = 11 dx[dy[ J I' dx'M W,,o{x[, y[; x',, y'^) x f{ Gt{x,\x'^)G^{y,\y'^) , (57) 

where W^,o(a;i, y'l; Xg, 1/2) is the Wigner function at t = and the Green's functions 
Gt{xj\x'j) are given in equation (fTH|) . The Wigner function W^fl{xi, yi',X2, y2) before the 
propagation is given by (remind that // = 0, a, and a is taken as real) 



exp 



" (27rai)(27ra2) 
with 

4 = ^(Aa)^^ • (59) 

A\ being as in equation (jlH)) . Since Wafi is Gaussian, the Wigner function Wa,t can be 
easily evaluated. One has 

„„„ J (xi+X2-iJ.e-^^y _ {yi+y2f _ {xi-X2-fJ.c-i^y _ {yi-y2? 
^^t* \ 4S^^ 4Ei 4S2^ 4S2_ 

= ^ (2.E'i)(2.E^_) 
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where 

The signals described by the Wigner functions of equation (pl|l correspond to entangled 
states if S2 < i 001 EH El , «-e. if 

e-2'' = l + Nx- ^Nx{Nx + 2) < e^* - (2M + l)(e^* - 1) . 

If no thermal noise occurs (M = 0), entanglement is present at any time, whereas for 
M 7^ the survival time is given by 



and the corresponding survival entanglement fraction /3s(A^, F, M, t) reads as follows 

The meaning of equation is that if the initial entanglement fraction is above 
threshold, i.e. [3 > /3s, the state remains not separable after propagation. Besides 
propagation noise, one should take into account detection efficiency at the heterodyne 
receiver. In this case the POVM of the detector is a Gaussian convolution of the ideal 
POVM 

n.(^) = A /d^/^ exp (-^^ j , (64) 



with 0"^ = (1 — ri)/iri. The Wigner function associated to the POVM (jMj) is given by 

W^[n,(ry)](xi,|/i;x2,2/2) = ^ "-^ "-^ (65) 

and the corresponding heterodyne distribution r^,r.j\f(^|«) by 

I2; — Q/ e~2^*p 

'^,;,r,A/(2;|a) = ^xl ^ X2 \ ' (66) 



with 



^lrM = 4S- + = (1 + 2Af)(l - e-^*) + e'^* + . (67) 

The error probability, as defined in the previous Sections, already taking into account 
that the optimal threshold is given by A = e~2'"*, reads as follows 



i?e(iV,r/,r,M) = ^ 



1 1 - Erf 


lae 2^ 


1 


2 A^^r.Af 









(68) 



and the optimal entanglement fraction /3opt(A^, J], P, M, t) = /^opt? which minimizes the 
error probability, is obtained after some algebra. One has 

2 jYo-2rt 

/^opt = , , (69) 

1 + A{N, r/, P, M, t) + B{N, r], P, M) v^l + C(A^, r/, P, M, t) 
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where 



A{N, r], r, M) = r/iVe-^*(2 - r] e"^*) + 4/(Ar, r/, T, M) , (70) 
B{N, r/, r, M) = 1 + r/iVe"^* + 2r/M(l - e"^*) , (71) 
C{N, r], r, M) = 2r]Ne-^\l - r] e"^*) + 2f{N, r/, T, M) , (72) 



with 



/(AT, r/, r, M) = 77M(1 - e"^*) [1 + r/A^e"^* + r/M(l - e"^*)] . (73) 

As a matter of fact, the optimal entanglement fraction depends on the noise parameters 
of the channel. In particular, it can be seen that (3opt decreases as F and M increase 
and 7] decreases. In other words, as the channel becomes more noisy, the entanglement 
becomes less useful. In figure El we report Popt for fixed detection efficiency and Ps as 
functions of the channel number of photons for different values of the other parameters. 

3.3. Direct vs heterodyne detection 

Heterodyne detection provides better performance than direct detection when the 
quantity 

D^Ar r 1 ^e(iV,^het,r,M) 

Be{N, r/ken, Vhet, T , M) = 1 - — — — — , (74) 

Ke{N, ?7ken,r, M) 

is positive, i.e. Re < K^. In equation (|7Hl ?7kcn and ?7het are the on/off and heterodyne 
detection efficiencies, respectively. In figure [7| we plot as a function of channel energy 
A^ for different values of the other parameters. Three regions of interest can be identified, 
since, as in the case of homodyne detection, for given values of r/ken, V^eti T and M we 
have two different thresholds for the energy channel, namely A^^*'j^°*'' (r/kcm ''7het; ^); 
j = 1,2, such that = 0. Heterodyne error probability is the smallest in a small 
region for A^ < N^^'^^^ and for large values of the channel energy (A^ > N^if), whereas 
in an intermediate interval of energy values ( A^th? '* < N < N^l^'^^^ ) the best results are 
obtained by direct detection (see table El). 

As for the homodyne detection, heterodyne receiver provides better results for either 
small or large values of the channel energy A^, whereas for the intermediate region, whose 
width depends on F, M and r] as in the case of A^, direct detection should be preferred. 

3.4. Heterodyne vs homodyne detection 

Since the channel energy intervals where heterodyne and homodyne detection should be 
preferred to direct one are quite similar, it is useful to introduce the function 

Ce(A/,r/,r/,l,M)-l ^^^^^^^^^^y (75) 

which is positive when < H^. As one can see in figure |Hl there exists a threshold 
Nth{t], r, M) on the channel energy A^ such that if A^ > A'th then Rc < He (see tablelH). 
Notice that the homodyne and heterodyne detection efficiencies have the same value. 
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In figure El we plot A'th for different physical situation. The threshold increases with 
increasing noise either in the propagation or in the detection stage. For a channel energy 
above the threshold the use of entanglement improves the communication performances. 

4. Conclusions 

In this paper we have analyzed binary communication in single-mode and entangled 
quantum noisy channels. We took into account different kind of noise that may occur, 
i.e. losses and thermal noise during propagation, non unit quantum efficiency of 
detectors during the measurement stage. 

As concern single mode communication, we found that, in presence of noise, 
homodyne detection is a more robust receiver when compared to direct detection. In 
particular, homodyne detection achieves a smaller error probability for either small or 
large values of the energy of the channel, whereas for intermediate values direct detection 
should be preferred. 

We then considered an entanglement based quantum channel build by amplitude 
modulated twin-beam and multiport homodyne detection. As for the homodyning, 
heterodyne detection should be preferred to direct detection for either small or large 
values of the energy of the channels. On the other hand the comparison between the 
performances of heterodyne and homodyne detection shows that there exists a threshold 
on the channel energy, above which the error probability using entangled channels and 
heterodyning is the smallest. The threshold depends on the amount of noise, and 
increases as imperfections in propagation and detection become more relevant. We 
summarized our results in tables HJ El El and El 

We conclude that entanglement is a useful resource to improve binary 
communication in presence of noise, especially in the large energy regime. 
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Channel energy 


Best detector (channel) 


N < 0.79 


Homodyne (single mode) 


0.79 <N< 4.46 


Direct (single mode) 


N > 4.46 


Heterodyne (entangled) 



Table 1. Comparison among direct, homodyne and heterodyne detection as a function 
of the channel energy when transmission and detection are ideal (see also figure EJ). 
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Channel energy 


Best detector (channel) 


N < iV,^7^ 


Homodyne (single mode) 


<N< 


Direct (single mode) 


N > <r^ 


Homodyne (single mode) 



Table 2. Comparison between direct and homodyne detection as a function of the 
channel energy in presence of noise during the transmission and detection stages. The 
thresholds A^th°™^(?7kcn, ^om, T, M), j — 1,2, define three diff'erent regimes (see also 
figure 131) • 



Channel energy 


Best detector (channel) 


N < 


Heterodyne (entangled) 


<N< 


Direct (single mode) 


N > 


Heterodyne (entangled) 



Table 3. Comparison between direct and heterodyne detection as a function of the 
channel energy in presence of noise during the transmission and detection stages. The 
thresholds N^^^*'\'r]]^cn,Vhct,^, M), j = 1,2, define three different regimes (see also 
figure [71). 



Channel energy 


Best detector (channel) 


iV<iVth 


Homodyne (single mode) 


iV>iVth 


Heterodyne (entangled) 



Table 4. Comparison between heterodyne and homodyne detection as a function of 
the channel energy in presence of noise during the transmission and detection stages. 
The threshold Nth{ri,r, M) defines two different regimes (see also figure |H1). 



in 



BS n 

' — -D 



ON/OFF 



reference 



Figure 1. Scheme of the receiver based on direct detection. The signal to be processed 
is mixed at a beam splitter of transmissivity r = cos^ with a given coherent reference 
state. Then, at the output, one of the two modes is ignored or absorbed, whereas the 
other one is revealed by on/ off photodetection. In the ideal case pin = Pf_i = 
/i = 0, a, and the the detection efficiency is = 1. In presence of noise, is given by 
equation l(TH|) and ij <1. 
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Figure 2. Scheme of the homodyne receiver. The signal to be processed is mixed at a 
balanced beam splitter with a highly excited coherent reference state, usually referred 
to as the local oscillator (LO). Then, the difference photocurrent is measured at the 
output. In the ideal case = p^ = = 0, a, and the the detection efficiency 

is ry = 1. In presence of noise, pin is given by equation H18f) and ry < 1. 




Figure 3. Plots of ylo(A^, rykcn, '7hom, T, Af), defined in equation (|4()|l . as a function of 
N for different values of Tt and M . The on/off and homodyne detection efficiencies 
are chosen to be the realistic values ryken = 0.95 and ?7hom = 0.85, respectively, and 
cos^ (f) — 0.99. The dashed line is with F, M — 0. In all the plots the solid lines 
correspond to (from bottom to top) M = 5 10^"^, 5 10~^, 0.1 and 0.2, respectively, while 
rt is: (a) r< = 10-2, (b) R = 5 lO'^, (c) Tt = 0.1, (d) Tt = 0.2. The insets refer to 
the region < TV < 1. When Ag > one has < K^. 



Binary communication in single-mode and entangled noisy channels 



19 



(a) 




(b) 




Type II 
non-collinear 





Figure 4. Schematic diagrams of possible implementations of entanglement based 
binary optical communication with TWB and multiport homodyne detection. In 
(a): TWB are produced by spontaneous downconversion in a type I nondegenerate 
optical parametric amplifier (NOPA). In this case the two modes have the same 
polarization. A displacement D{a) is applied to mode a and then the signal is 
transmitted toward a multiport homodyne receiver. In (b): TWB are produced by 
spontaneous downconversion in a noncoUinear type II NOPA. In this case the two TWB 
modes have different wave-vector and orthogonal polarizations. A displacement D{a) 
is applied to mode a, while a quarter wave plate (QWP) rotates the polarization of the 
mode h so that they become same polarized. After the propagation, the two modes 
are detected as before. In (c): TWB are produced by spontaneous downconversion 
in a coUinear type II NOPA. In this case the two TWB modes have the same wave- 
vector and orthogonal polarizations. A first polarizer beam splitter (PBS) reflects the 
mode a, which is displaced by -D(a), and transmits mode &; the two modes are then 
recombined in a second PBS and transmitted. Before detection another PBS reflects 
mode a toward a QWP that changes its polarization to the same of mode h. During 
propagation modes are subjected to losses (at rate P) and thermal noise (parameter 
M). The detector efficiency is 77. 




Figure 5. Linear-log plot of the error probabilities of direct (Kc, dotted line), 
homodyne {Ho, dot-dashed line) and heterodyne {Re, solid line) detection in the ideal 
case (absence of noise in propagation and detection) as functions of the channel energy 
N. We chose cos 0^ = 0.99 for direct detection. As the channel energy increases, the 
heterodyne error probability decreases more rapidly than the single-mode ones. In 
the smaller pictures: magnification of the plot in the regions (a): < iV < 1 and 
(b): < iV < 5. For a channel energy TV < 0.79 homodyne detection gives the best 
performances, whereas as the channel energy increases the best results are obtained 
by the direct (0.79 < < 4.46) and heterodyne detection {N > 4.46). 
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Figure 6. Plots of the optimal entanglement fraction /?opt (solid lines) given in 
equation H69() and of the survival entanglement fraction /3s (dashed lines) given in 
equation 1)63(1 as a function of the channel energy N for different values of F, M and 
rj. The plots (a2) and (b2) are magnifications of the regions < iV < 0.2 of plots (al) 
and (bl), respectively. In all the plots we put (from top to bottom for solid lines, from 
bottom to top for dashed lines) M = 5 10"'^, 5 10~^, 0.1 and 0.2, respectively, whereas: 
(al) 7/ = 0.9, Tt = 10"^ (bl) r] = 0.9, Ft = 210"^ As the channel becomes more 
noisy, entanglement becomes less useful. Notice that for this choice of the parameters, 
when N > 0.056 (a2) and N > 0.142 (b2) the state arriving at the receiver is always 
non separable. 
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Figure 7. Plots of _Bc(^, '7ken, ^ot, 1^7 ^^), defined in equation 1(74(1 . as a function of 
N for different values of Ft and M. The on/off and heterodyne detection efficiencies 
are chosen as 7?ken = 0.95 and ?7het = 0.85, respectively, and cos^ (j) = 0.99. The dashed 
line is with F, M = 0. In all the plots the solid lines correspond to (from bottom 
to top) M 5 10"^, 5 10^2, 0.1 and 0.2, respectively, whereas Ft is: (a) Ft 10^^, (b) 
Ft = 5 10-^ (c) Ft = 0.1, (d) Ft = 0.2. The insets refer to the region < iV < 0.5. 
When i?e > one has Re < Kg. 
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Figure 8. Plots of Ce{N,ri,T, M), defined in equation H75|l . as a function of N for 
different values of Tt and M. The heterodyne and homodync detection efficiencies are 
chosen as 77=0.85. The dashed line is Co with F, AI = 0. In all the plots the solid 
lines correspond to (from top to bottom) M = 5 10^"^, 5 10^^, 0.1 and 0.2, respectively, 
whereas Tt is: (a) Tt ^ 10"^, (b) Tt ^ 5 10'^, (c) Tt = 0.1, (d) Tt = 0.2. When Cc > 
one has Rg < Hg. 
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Figure 9. A^th as a function of Tt for different values of the detection efficiency 77 and 
the average number of thermal photons M. When N > N^^ the error probability using 
entangled channels and heterodyne detection is less than the one obtained with a single- 
mode homodyne detection (i?c < He). The dashed line is A^th for 77 = 1 and M = 0, 
while the solid lines represent the threshold for (from right to left) M = 0, 5 10^^, 0.1 
and 0.2, respectively, whereas: (a) = 1, (b) 7; = 0.9, (c) 77 = 0.8 and (d) 77 = 0.7. 



